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ABSTRACT 

We present an analysis of the neutron star High Mass X-ray Binary (HMXB) 4U 1909+07 mainly based on Suzaku data. We extend 
the pulse period evolution, which behaves in a random-walk like manner, indicative of direct wind accretion. Studying the spectral 
properties of 4U 1909+07 between 0.5 to 90keV we find that a power-law with an exponential cutoff can describe the data well, when 
additionally allowing for a blackbody or a partially covering absorber at low energies. 

We find no evidence for a cyclotron resonant scattering feature (CRSF), a feature seen in many other neutron star HMXBs sources. 
By performing pulse phase resolved spectroscopy we investigate the origin of the strong energy dependence of the pulse profile, 
which evolves from a broad two-peak profile at low energies to a profile with a single, narrow peak at energies above 20keV. Our 
data show that it is very likely that a higher folding energy in the high energy peak is responsible for this behavior. This in turn leads 
to the assumption that we observe the two magnetic poles and their respective accretion columns at different phases, and that these 
accretions column have slightly different physical conditions. 
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1. Introduction 

High mass X-ray binaries (HMXBs), in which the compact ob- 
ject is a rotating neutron star, provide a unique way to investi- 
gate the flow of matter into a deep gravitational potential. By 
analyzing the evolution of the pulse period and its correlation to 
luminosity, the presence or absence of a stable accretion disk can 
be probed (Ghosh & Lamb, 1979). By performing pulse phase 
resolved spectroscopy, different areas of the X-ray producing re- 
gion can be analyzed and thereby the interaction between ra- 
diation and matter in very strong magnetic fields can be stud- 
ied. 4U 1909+07 is an ideal source for these investigations, as 
it is a slowly rotating, persistent HMXB. Discovered by Uhuru 
(Giacconi et al., 1974), its 4.4 d orbital period was determined in 
RXTE/ASM data by Wen et al. (2000). Pulsations with a pulse 
period of «605 s were later found in RXTE/PCA data by Levine 
et al. (2004). These authors also found that the intrinsic absorp- 
tion Afa is strongly variable over the orbit, peaking around orbital 
phase 0. These measurements can be explained by a smooth stel- 
lar wind and an inclination of the system between 38° and 68° 
(Levine et al., 2004). Fiirst et al. (2011a, hereafter Fll) investi- 
gated the long-term evolution of the pulse period between 2003 
and 2009 using INTEGRAL flSGRl lightcurves. We found that 
the evolution is consistent with a random walk, a typical behav- 
ior for a wind accreting source without a stable accretion disk, 
similar to Vela X-l (see, e.g., de Kool & Anzer, 1993). 

The spectrum of 4U 1909+07 can be described using typical 
phenomenological models often applied to HMXB neutron star 
sources, such as a power-law with an exponential cutoff (Levine 
et al., 2004, Fll). By analyzing phase-resolved RXTE data we 


found that a blackbody component with a temperature around 
1.7 keV is also needed to describe the soft excess in the spec- 
trum. Many neutron star sources show a cyclotron resonant scat- 
tering feature (CRSF, see, e.g., Schdnherr et al., 2007) at en- 
ergies between 10 and 50 keV. Neither in RXTE/HEXTE nor 
in INTEGRAL fIS GRI data is evidence for such a feature seen; 
however, the signal-to-noise ratio (S/N) in these data is rather 
low. To investigate the behavior of the soft excess in more de- 
tail and to perform a rigorous search for a CRSF, we obtained a 
20 ks Suzaku observation, the analysis of which is presented in 
this article. In Sect. 2 we give a short overview of the data used 
and the reduction pipeline. In Sect. 3 we discuss the evolution of 
the pulse period, while in Sect. 4 we analyze the X-ray spectrum. 
We perform phase averaged (Sect. 4. 1) as well as phase resolved 
spectroscopy (Sect. 4.2). In Sect. 5 we summarize and discuss 
our results. 


2. Observations & data reduction 

This paper is mainly focused on a 20 ks Suzaku observation per- 
formed on November 2, 2010. Data were taken in the orbital 
phase range 0.472-0.627, i.e., outside the phases of increased 
absorption. We use data from the X-ray Imaging Spectrometers 
(XISs, Koyama et al., 2007) 0, 1, and 3 in the energy range 0.5- 
lOkeV and from the High X-ray Detector’s (HXD, Takahashi 
el al., 2007) PIN instrument from lOkeV to 90keV. We did not 
use data from the HXD/GSO due to source confusion with the 
nearby, bright black-hole binary GRS 1915+105. As the hard 
X-ray spectrum of 4U 1909+07 falls off exponentially, a sig- 
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Fig. 1. Lightcurves taken with XIS 0 (red), X1S 1 (orange), and XIS 3 
(green) in the 0.5-10 keV energy band and PIN (blue) between 10 to 
29keV, rebinned to 160 s. 


nificant detection in GSO would be rather unlikely. The data 
were reduced using the standard pipeline from HEASOFT v6. 1 1 
and the HXD calibration from September 2011, including dead- 
time correction for PIN. Data analysis was performed using 
the Interactive Spectral Interpretation System (ISIS, Houck & 
Denicola, 2000), version 1.6.2. Wc extracted lightcurves from 
XIS with 16 s time resolution in the energy bands 0.5-5 keV, 
5-10 keV, and 0.5-10 keV, as well as from PIN with 40 s res- 
olution in the energy bands 10-20keV, 20-40 keV, and 40- 
90 keV. Figure 1 shows the lightcurves from all instruments, re- 
binned to 160 s for better visibility. Some flaring is evident in the 
lightcurve, however, no significant changes in hardness could be 
found when comparing the different energy bands. For compar- 
ison we also extracted data from the ISGRI detector (Lebrun 
et al„ 2003) aboard INTEGRAL (Winkler et al., 2003), between 
MJD 55107.13 to 55148.82 and MJD 55503.74 to 55518.02 
with exposure times of 235ks and 396 ks, respectively. Using 
the standard OSA 9.0 we extracted lightcurves between 20 to 
40 keV with 20 s time resolution. All lightcurves, Suzaku and 
INTEGRAL, were transferred to the barycenter of the solar sys- 
tem and corrected for the orbital motion of 4U 1909+07, using 
the ephemeris by Levine et al. (2004). 


3. Timing 

The pulse period evolution of 4U 1909+07 was published by 
FI 1 up to spring 2009. Figure 2 extends this time series up to 
fall 2010, using all publicly available INTEGRAL data and the 
Suzaku observation. INTEGRAL also observed 4U 1909+07 
in spring 2010, but the rough sampling does not allow for 
a reliable pulse period measurement. To obtain the pulse pe- 
riods we performed epoch folding (Leahy el al., 1983) us- 
ing XIS lightcurves between 0.5-10 keV with 16 s time resolu- 
tion, PIN lightcurves between 10-90 keV with 40 s time reso- 
lution, and INTEGRAL / ISGRI lightcurves between 20-40 keV 
with 20 s time resolution. Uncertainties on the pulse periods 
of INTEGRAL were obtained by fitting a Gaussian distribu- 
tion to the x 1 distribution of the individual epoch folding re- 
sults, as described by Fll. The epoch folding distribution of 
Suzaku/XIS as well as Suzaku/PIN does not show a well enough 
defined peak for fitting with a Gaussian function. We therefore 
estimated the uncertainties on the period measurement using a 
Monte Carlo approach as described by Davies (1990), where 
we simulated 1000 lightcurves with the same binning as the 
original measurement, scattered around the measured pulse pro- 
file. The width of the distribution of resulting pulse periods in 
these lightcurves is a good indicator of the uncertainty of the 
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Fig. 2. Pulse period evolution between 2003-2011. Black data were 
already published by Fll. More recent INTEGRAL /ISGRI data are 
shown in red, Suzaku/PIN in blue. Suzaku/XIS is not shown to ensure 
readability. 


pulse period. This approach leads to Pxis = 604. 10 ± 0.20 s and 
Ppin = 603.87 ± 0.24 s (all uncertainties are given at the 90% 
level). Due to the longer exposure time, the almost simultane- 
ous INTEGRAL/ISGRT data provide a much better constrained 
pulse period measurement, with Pisgri = 603. 86± 0.05 s. As this 
period is consistent with the Suzaku measurements, we will use 
it throughout this paper and for the analysis of the Suzaku data. 

The pulse profiles of Suzaku/PIN and INTEGRAL/ISGRI 
in the 20-40 keV range agree, showing a short single peak fol- 
lowed by a long dim phase (Fig. 3d). This profile is in agree- 
ment with RXTE/HEXTE (see, e.g., Levine et al., 2004, FI 1). 
As described by Fll a strong change in the pulse profile occurs 
between 8 and 20 keV, where it changes from a two-peaked pro- 
file to a single peaked profile. This energy dependent behavior is 
also seen in other neutron star X-ray binaries, like 4U 01 15+634 
(Muller et al., 2010) and 1A 1118-61 (Suchy et al., 2011). The 
energy range of the transition is only marginally covered by the 
instruments aboard Suzaku, as it is just in the overlapping re- 
gion between XIS and PIN. The XIS pulse profile in different 
energy bands is shown in Fig. 3a and Fig. 3b and is very similar 
to RXTE/PCA profiles as presented by Levine et al. (2004) and 
Fll. 


4. Spectroscopy 

4.1. Phase-averaged spectroscopy 

As discussed by Fll, the phase averaged spectrum of 
4U 1909+07 can be very well described with typical phe- 
nomenological models often used to describe the spectra of neu- 
tron stars in HMXBs plus a Gaussian-shaped fluorescent Fe Kor 
line. The XSPEC models cutoffpl, highecut (White et al., 
1983), and comptt (Hua & Titarchuk, 1995) all gave similar 
statistically acceptable fits. Therefore, we applied these mod- 
els to the Suzaku data and the results are shown in Tables 1, 2, 
and 3, respectively. The data are shown together with the best-fit 
cutoffpl model in Fig. 4. 

The flux of the Suzaku data (XIS normalization 1 ) is compa- 
rable to the RXTE data (PCA normalization) presented bv FI 1, 
= 0.324!°™ keV S -' cm- 2 and = 

0.330 ± 0.007 keV s" 1 cm 2 , respectively. The Fe K a line is nar- 
row and its width could not be resolved in the XIS data. We 

1 PIN cross-calibration factor was within 2% 
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Fig. 3. Energy resolved pulse profiles as measured with SuzakuptlS (a) 
and b)) and PIN ( c ) to e)). Marked in color and labeled A-E are the 
phase bins used for phase-resolved spectroscopy (see text for details). 
In panel c), INTEGRAL/1SGRI data is superimposed in red, scaled ac- 
cording the right-hand y-axis. 


therefore fixed it to 10“ 6 keV in all fits. It was shown by FI 1 
that an additional blackbody component is necessary to describe 
a soft excess. Adding this component leads also to a clear im- 
provement in the fit to the Suzaku data, as can be seen in the 
second column of Tables 1, 2, and 3. However, some continuum 
parameters become only very weakly constrained, like Ef oW in 
the cutoffpl and highecut model or the seed photon tem- 
perature 7o in the comptt model. Additionally the 20-40 keV 
flux drops to unrealistically low values 2 and the best fit black 
body temperature is very different between the cutoffpl and 
highecut model, rendering a physical interpretation impossi- 
ble. To describe the soft energy part of the spectrum without 
a black body component, we included a partial coverer in the 
model with two different absorption columns Nm and 7 Vh i2 , with 
their relative influence determined by a covering fraction CF: 

(CF X Afa.i + (1 - CF) x Ar Hi2 ) X (cont + Fe Ka) (1) 

This model leads to equally large improvements in terms of* 2 , 
as seen in the third column of Tables 1, 2, and 3. Although 
the secondary absorption column is also only weakly con- 
strained, the continuum parameters remain very well constrained 
and at physically reasonable values. 

4 . 2 . Phase-resolved spectroscopy 

To study the variance of the spectrum with the viewing angle 
onto the neutron star, i.e., the pulse phase, we divided the data 

- compensated by a very high PIN cross-calibration factor of » 1.3 


Table 1. Phase averaged spectral parameters for the cutoffpl con- 
tinuum model. In the second column the best-fit parameters with an 
additional blackbody component (bbody) are shown, in the third the 
parameters when including a partial coverer (PC). 


Model 

cutoffpl 

cutoffpl 

cutoffpl 

parameter 


+bbody 

+PC 

^3-!0tev“ 

9.33 ±0.07 

9.1*7 ± 0.07 

9.29 ±0.07 

?4M0ktV° 

8.6 ±0.7 

97* \ 
*- 7 -0.4 

8.7 ±0.7 

Nh-1* 

JV h -2‘ 

CF 

10.71 ±0.30 

10.84:1 

0,25 ±0.06 

r 

0.93 ±0.06 

1.09^3 

22 ™ 

fi 'IRR+O - 016 
°- J88 -0.005 

1.20 ±0.09 

Em DkeV] 
Et* [keV] 

15.4 +t7 

6.388!^ 

19.3+™ 

6.391lgSg 

A Fe' 

1.98 ±0.18 

1.84 ±0.18 

1.94 ± 0.20 

MbbtkeV] 

- 

2.03^g 

- 

Abb 

- 

(0.98:§j»)xl0- 3 


X 2 /iof 

544.78/414 

506.74/412 

497.76/412 


Notes. (a) in 10 “ 2 keV s -1 cm “ 2 W in units of 10 22 H- atoms cm -2 
(c) Covering Fraction in I0 " 4 phs 1 cm -2 . 


Table 2. Same as Tab. 1, but for the highecut model. 


Model 

highecut 

highecut 

highecut 

parameter 


+bbody 

+PC 

T^lOkev" 

9.28 ± 0.07 

9.27 ±0.07 

9.28 ± 0.07 

7*20-40 keV* 

Afe-1* 

8.1 ±0.7 
10.3 ± 0.4 

14.551 

8-6$ 

9 o*i*ro 

N H -2 b 

CF 

— 


i°- 4 ^ 

0.15*®X 

r 

Efoid [keV] 

1 - 12:88 

i8.o:{| 

1-3*58 

24™ 
70 ij 

, .aiftlS 
6-390% 

1-963:5. 

Ecui [keV] 

6.0 ±04 

Er t [ keV] 
A*' 

6390t8g 

1.82 ±0.19 

6.3839^ 

1 qo+Oj 

of 4 xn 

U J ;-aio 

&»>[ keV] 



Abb 

- 

(2.3:jf) x 10~ 3 

- 

A^/dof 

499.20/413 

494.65/411 

492.03/411 

Notes. (fl) in 

10 _2 keV s" 1 

cm 2 (b) in units 

of lO^H-atomscr 

(c) Covering Fraction in 10 4 ph s -1 cm -2 . 


Table 3. Same as Tab. 1, but for the compTT model. 

Model 

compTT 

compTT 

compTT 

parameter 

nr a 

0 nn ^ n nn 

+bbody 

no-r T n nn 

+PC 

rv . A rv-> 


^20-40 keV fl 

Nh- 1 * 

N k - 2 * 

CF 

7b[keV] 
kT [keV] 
r 

[keV] 

Afc" 

*Tbb[keV] 

Abb 

*Vdof 


8.6 ± 0.8 
6.3 ± 0.4 


1.18 ±0.06 
n a+O.7 


67 tl 

io.3:f 0 


O 04 + 0.15 


9.4^ 

7 “* 

7 , 1 - 0.6 

0.28 ± 0.06 


OS 

4.3 ±0.4 

6-388% 
1.83 ±0.18 

/O -0.5 

4.1 ±0.4 

6.388!°™ 
1.87 ±0.17 

6.387!$“ 
1 qq+ozst 
1 V *- 0.20 

— 

2.13 ±0.13 


- 

(1.41 ± 0.20) xlO - 3 

- 

551.26/413 

501.33/411 

506.79/411 


Notes. W in 10 2 keV s ’" 1 cm " 2 {b) in units of 10 22 H-atomscm " 2 
(c) Covering Fraction [d) in 10 " 4 ph s 1 cm -2 
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Energy [keV] 

Fig. 4. a) XISO (red), XIS1 (green), XIS3 (orange), and PIN (blue) spec- 
trum together with the best fit cutoffpl-model including partial cov- 
ering. In magenta the strongly absorbed part is shown, in dark blue the 
less absorbed one. The Fe K a line is shown in brown, b) residuals to 
the cutoffpl model without partial coverer. c) residuals to the best-fit 
model. 


in 5 phase bins, as indicated in Fig. 3. The bins were chosen 
to cover the prominent features, e.g., the deep minimum be- 
tween the pulses and the hard primary pulse (bins C and D, 
respectively) as well as to retain a S/N sufficient to constrain 
the individual spectral parameters well. As shown by Fll, the 
cutoffpl model delivers the best description of the phase- 
resolved spectra, and in order to be able to compare our results 
directly with the previous work, we applied that model to the 
phase-resolved data. Neither a blackbody component nor a par- 
tial coverer could be significantly detected in these spectra due to 
the lower S/N, therefore we only use a single absorber, the con- 
tinuum, and a Gaussian-shaped fluorescent iron line at 6.4 keV. 

Detailed analysis of the phase averaged spectrum shows that 
there are strong systematic correlations between Ah, T, and £f 0 id. 
These correlations lead to large uncertainties in the phase re- 
solved spectra in all parameters. To reduce these uncertainties 
we fixed each of the three parameters in turn, and compared the 
results with the energy resolved pulse profiles. All fits lead to 
very similar ;^ d values, so that a distinction based on the quality 
of the fit is not possible. However, when comparing the fluxes 
in the 20-40 keV energy band modeled with a fixed folding en- 
ergy to the pulse profile in that energy range, a clear discrep- 
ancy is evident, i.e., the bright main peak is not correctly re- 
produced. Additionally, Fll found indications that the phase re- 
solved RXTE spectra are best described using a variable E fold- 
We therefore exclude the model with a fixed folding energy from 
the following discussion. 

As can be seen in Fig. 5f, the folding energy is in all mod- 
els very high during phase bins C and D, the minimum and the 
primary peak. The folding energy seems therefore responsible 
for the spectral hardening during the primary peak, and not, e.g., 
strong changes in the photon index T. The results also show that 
the minimum has a similar spectrum as the primary peak. In the 
models with variable T, £foid stays high during phase bin E, but 
at the same the time T is rising, leading to an overall softer spec- 
trum. With T frozen, £f 0 id drops to values similar to those in 
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Fig. 5. Results of the phase resolved spectroscopy. In black the fit with- 
out any fixed parameters is shown, for the blue data points the A H was 
frozen, for the red the photon index T. In green the results of the phase 
averaged fit are shown, a) XISO pulse-profile in the 0.5-10 keV range, 
b) flux between 3-10keV in keV s _1 cm" 2 , c) flux between 20-40 keV 
in keV s 1 cm -2 , d) absorption column density, e) photon index,/) fold- 
ing energy, g) flux of the iron line in phs~ l cm -2 . 


phase bins A and B. Differences between these two descriptions 
would only be visible above ~70keV, where the source flux has 
dropped below the detection limit of PIN. 

Variations of Ah are visible in Fig. 5d, but seem uncorre- 
lated with the observed hard primary peak. Ah drops to low val- 
ues during the minimum in phase bin C, but the overall variation 
is only marginally significant. It is interesting to note that the 
flux of Fe Ka is also variable with phase (Fig. 5g), however it is 
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not correlated with the overall soft X-ray flux nor the Afa varia- 
tions. There rather seems to be a slight shift with the Fe Ka flux 
reacting on the 3-10 keV flux with approximately 120 s delay. 
Using the model in which all parameters were allowed to vary 
aid shifting the Fe Ka flux by one phase bin, the correlation 
to the soft continuum increases from 0.04 to 0.70 according to 
Pearson’s correlation coefficient. The energy of the Fe Kcr line 
(not shown) is not varying significantly with pulse phase. 

5. Summary & Outlook 

We have presented a timing and spectral analysis of a 20 ks 
Suzaku observation of 4U 1909+07. We extended the pulse pe- 
riod evolution as presented by Fll by two data points, the 
most recent of which, P = 604.86 s, is consistent between 
INTEGRAL and Suzaku. To describe the spectrum we applied 
different phenomenological models and find statistically accept- 
able fits when including a fluorescent Fe Ka line at 6.4 keV and 
either a blackbody component or a partially covering absorber. 
The continuum parameters in the models with a blackbody com- 
ponent are very similar to the ones found by FI 1 in RXTE- 
data. However, they are highly unconstrained in the Suzaku- 
data, while using a partial coverer still results in a good descrip- 
tion of the continuum with small uncertainties. 

A partial coverer can be used as an approximation to a com- 
plex structured stellar wind (“clumpy” wind), as seen in many 
HMXBs as well as isolated O and B stars (see, e.g., Oskinova 
et al., 2007; Fiirst et al., 2011b, among many others). Within 
these winds, dense structures with column densities of the order 
of 10 24 cm -2 can be reached, as seen in many other sources (e g., 
GX 301-2 or IGR J16318-4848; Fiirst et al., 2011b; Barrage 
et al., 2009). The absorption column of the partial coverer in 
4U 1909+07 is measured to be also of that order, with values 
between (O.S-lJxlO 24 cm -2 , as seen in Tables 1, 2, and 3. 

We also applied the partial covering model to the older 
RXTE data, which have about twice the column density as the 
Suzaku data, but did not find an acceptable fit. As the RXTE 
data do not cover the very soft energies, a small contribution of 
the partial coverer can easily be missed in the data. If the partial 
coverer is due to the stellar wind, it is also highly variable, as 
s:ructures will move with the wind with typical outflow termi- 
nal velocities > 1000 km s -1 (Prinja et al., 1990). Therefore, the 
RXTE observation might have just been performed during a less 
obstructed view onto the X-ray source. Nonetheless, these data 
c’early require a blackbody component to describe the data. 

Assuming that the model including the partial coverer repre- 
sents a good description of the physical conditions in and around 
the X-ray source, it means that the blackbody disappeared, i.e., 
that it is highly variable. The blackbody component originates 
from the thermal mound on the neutron star or the accretion col- 
umn and its properties can vary with changes in the accretion 
rate and the physical conditions in the accretion column close 
to the neutron star. The X-ray fluxes of the Suzaku and RXTE 
data are, however, similar, i.e., it is likely that the accretion rate 
and thereby the temperature in the column are also very simi- 
lar. Changes in the geometry of the accretion column could ac- 
count for small differences, like the ones seen when comparing 
the RXTE model to the best fit Suzaku model which includes a 
blackbody component. In Suzaku the temperature is measured to 
be around 2 keV, compared to « 1.5 keV in the RXTE data. That 
the blackbody component becomes completely invisible due to 
changes in the geometry of the accretion column seems unlikely. 
It is also possible that the blackbody is a spurious effect in the 
RXTE data. For example, the galactic ridge emission, which can 


strongly influence RXTEfPCA data due to the large field of view 
of the collimator, can also often be described by a soft thermal 
component (see, e.g., Yamauchi et al., 2009). The galactic ridge 
component is suppressed in Suzaku/XIS due to its imaging capa- 
bilities. To study the variability of the blackbody and the partial 
coverer, a more extensive monitoring in the soft X-rays is neces- 
sary. 

To investigate the behavior of the spectrum under different 
viewing angles, we performed phase resolved spectroscopy in 5 
phase bins. As already seen in the energy resolved pulse profiles 
(Fig. 3), the primary peak (labeled D) is distinctly harder than the 
broad secondary maximum (bins A, B, and E). We found that 
this is likely due to a reduced folding energy in the secondary 
maximum. This change can be explained by assuming that ac- 
cretion happens on both poles of a magnetic dipole field and 
that each pulse maximum corresponds to one of the poles. Small 
differences in magnetic field strengths will lead to different ge- 
ometries and thereby temperatures in the two accretion columns. 
The hotter the electron plasma the more effectively it can scatter 
photons to high energies via the inverse Compton effect, leading 
to an overall harder spectrum. This picture might be over sim- 
plified, as relativistic light bending effects and the geometry of 
the accretion column need to be taken into account, but such an 
analysis is beyond the scope of this article. There are, however, 
new models under development to describe the X-ray producing 
region more realistically (e.g., Becker et al., 2012), which we 
will apply in a future publication. 

We searched for CRSFs in the phase averaged as well as in 
the the phase resolved spectrum, but did not find any evidence 
for such a feature in the energy range between 10 and 70keV. 
Assuming typical values of Ecrsf = 25keV and ocrsf = 
10 keV the maximal depth for such a line would be tcrsf ^ 0-23, 
assuming a line with a Gaussian optical depth. This is roughly a 
factor of two less than the depths of the lines seen for example 
in A 0535+269 (Caballero et al., 2007) or 4U 1907+09 (Rivers 
et al., 2010). As 4U 1909+07 is clearly a young neutron star, 
as proven through the strong pulsations and the early type stel- 
lar companion, it also possesses a strong magnetic field, which 
in theory must lead to the formation of CRSFs. As theoretical 
simulations show, however, it is possible that these features (es- 
pecially die fundamental line) can be filled up again by photon 
spawning and are therefore not detectable in the observed X-ray 
spectrum (Schdnherr et al., 2007). Additionally, the geometry of 
the magnetic field as well as the viewing angle onto the neutron 
star both can have strong effects regarding line shape and depth. 
Work is still in progress in understanding and quantizing these 
effects. The stringent upper limit on the CRSF depth is therefore 
very important to better understand the physics and magnetic 
fields of neutron stars. 
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